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ABSTRACT 

We summarize observations with the Spitzer Infrared Spectrograph (IRS) of 
571 starbursts (strong PAH emission features), 128 obscured AGN (strong sil- 
icate absorption), and 39 unobscured AGN (silicate emission). Sources range 
in luminosity from 10* to 10^'^ L0 and continuously in redshift for < z < 3. 



O 
(N 

>>■ 

cn 

O 

u 

Qh' The most luminous starbursts and AGN evolve as (1+z)^'^ to z ~ 2.5; no clear 

O ■ evidence is found that this evolution ceases beyond z = 2.5. Dust obscuration 

^ • in starbursts is determined by comparing PAH luminosity with ultraviolet lumi- 

, ^, ' nosity and indicates severe obscuration in most starbursts, even those selected 

in the ultraviolet; the median ratio (intrinsic ultraviolet /observed ultraviolet) is 
^ ■ ~ 50 for infrared selected starbursts and ~ 8 for ultraviolet selected starbursts. 

^ ■ Obscuration increases with bolometric luminosity, but starbursts which appear 

C<S ■ most luminous in the ultraviolet are those with the least obscuration. This result 

indicates that extinction corrections are significantly underestimated for ultravi- 
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o, 

O . deduced only from rest frame ultraviolet observations 



Q ■ olet selected sources, suggesting that galaxies at z > 2 are more luminous than 
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Vh ' 1. Introduction 



At the same time that the Spitzer mission was being conceived as NASA's final Great 
Observatory (originally as SIRTF), results from the Infrared Astronomical Satellite (IRAS) 
demonstrated that the most luminous galaxies in the universe are dusty and heavily ob- 
scured, with most of their luminosity arising from optical and ultraviolet luminosity which 



is reradiated by the absorbing dust ( ISoifer. Neugebauer. fc Houckl 119871 ). A fundamental 



motive of Spitzer was to track the evolution of these luminous galaxies through the history 
of the Universe. In what follows, we review the success of Spitzer in achieving this goal. 

Our conceptual sketch of such an ultraluminous infrared galaxy (ULIRG) is shown in 
Figure 1 and described in the caption. Undoubtedly, starbursts are present to some extent 
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in all such sources. The starburst ULIRGs are the most fertile star forming regions in the 
Universe. However, the observed characteristics are determined both by the luminosity of 
the starburst and the luminosity of dust heated by a central active galactic nucleus (AGN). 
Here, we separate the two phenomena to track separately the evolution of starbursts and 
AGN. 

The plots and results that we show below derive only from IRS low resolution obser- 
vatio ns. More detailed diagnostics can be made for those sources with high resolution data 
(e.g. iFarrah et al.l 120071 : iBernard-Salas et al.l 120091 ). but the heavy extinction in the most 
obscured sources can hide even infrared emission lines arising near the central AGN. For ex- 
ample, the presence of only blueward asymmetries on infrared forbidden lines from heavily 
obscured ULIIRG AGN in dicates that the obs curing dust is thick, even at the scale of the 
narrow line region of AG N (ISpoon &: Holtl 20091 ). In general, the new IRS results confirm the 



initial demonstration by iGenzel et al. 



( 19981 ) that strong emission from polyclyclic aromatic 
hydrocarbons (PAH) is associated with starbursts, but not with AGN. The photodissociation 
region at the interface between the surrounding molecular cloud and the HII region of a star- 
burst provides the a ppropriate physical conditions for strong PAH excitation and emission 
(jPeeters et al.ll2004J ). The hundreds of low resolution spectra accumulated by the IRS give 
an empirical result that "pure" starbursts - sources with no evidence in any spectral region 
of an AGN - always have equivalent width for the 6.2 ^m PAH feature, EW(6.2 /im), to 
be EW(6.2 um)> 4 /xm JSrandl et al.lbood : iHouck et al.lboOTi : IWeedman fc Houck|[20093 : 
Sargsyan et al.ll2010l ). 



Conversely to starbursts, there are large numbers of sources with IRS spectra which 
are "pure" AGN, having no evidence in any spectral region of an observationally significant 
starburst; these have EW(6.2 /im)< 0.1 /im. We have used these simple observat ional criteria 



to der ive t he samples of starbursts and AGN presented below, as explained in iHouck et al. 



( 20071 ) and lWeedman fc Houckl ( 20081 . l2009a| ]bl). Our primary motivation for using the simple 
observational measurement of peak flux densities for the PAH 7.7 fim feature in spectra of 
"pure" starbursts or the 7.9 /im spectral peak for "pure" obscured AGN is to have reliable 
and self-consistent data for the large numbers of sources discovered by Spitzer at z > 2, for 
which S/N is poor, and deconvolution of the total flux within spectral features is uncertain. 

To measure total infrared luminosities, we have used as an initial calibration a limited 
number of ULIR Gs or starbursts hav i ng IR AS fluxes and IRS spect ra, so that Lhoi can 
be estim ated as in ISanders fc Mirabel I (119961 ). These calibrations, from lWeedman fc Houck 
( l2009bl ). are given in the captions of Figures 3 and 4. This estimate of L^o; arises only from the 
IRAS wavelengths, with an assumed correction for colder dust components, so our calibration 
will be improved as more far infrared and submillimeter observations accumulate. We use 



-3- 









k 



.005 loGi ^ 



Fig. 1. — Sketch of the assumed model for a luminous, dusty galaxy. Starbursts can be found 
within the large torus containing molecular clouds and dust; each starburst immersed in a 
molecular cloud is surrounded by a photodissociation region (PDR) which produces strong 
PAH emission features. Any active galactic nucleus (AGN) heats dust close to the nucleus. 
(Note logarithmic size scale at bottom.) The starbursts will appear the same from any 
viewing direction, but the appearance of the AGN depends on whether it is viewed face on 
or through the torus. If face on, the hottest dust is seen directly and silicate emission is seen. 
If seen through the torus, the dust continuum from the hot dust is absorbed within cooler 
dust and produces sihcate absorption. The final observed spectrum is some combination 
of PAH, silicate absorption, or silicate emission depending on the relative luminosities of 
starbursts and AGN, and on the viewing direction for the AGN. 



this L hni to cahbrate the star formation rate (SFR) according to the precepts of iKennicutt 



( Il998l ). as given i n the captioii to Fig ure 3. All sources shown below were observed with the 
IRS, described in lHouck et al.l (|2004[ ). using only the low resolution modules, giving spectral 
coverage from ~5 /xm to ~35 /im. For all sp ectra which we an alyzed, spectral extractions 
were done with the SMART analysis package (JHigdon et al.ll2004j ). usually with the improved 
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Fig. 2. — Representative rest frame spectra observed for the most luminous dusty galaxies 
at z = 2.5. The top spectrum represents the brightest type 1 AGN and is featureless at these 
rest wavelengths, but silicate emission begins to increase at the longest wavelengths toward a 
peak at ~ 10 ^m. Such sources are optically bright and have known redshifts, so IRS spectral 
measures of the dust continuum are made at 7.9 //m. The center spectrum is a heavily 
obscured AGN with an apparent spectral peak at 7.9 //m arising because of deep silicate 
absorption at longer wavelengths and other absorption features at shorter wavelengths. The 
dust continuum is measured at this 7.9 fim peak. The lowest spectrum is a starburst with 
strong PAH features. Starburst luminosity is measured by the peak of the strong 7.7 /im PAH 
feature. (The illustrated spectra have higher S/N than individual sources actually observed 
at such redshifts; spectra which are shown derive from brighter, lower redshift examples.) 



"optimal extraction" procedure in iLebouteiller et al.l ( 2010[ )lH 



^http://isc. astro. cornell.edu/IRS/Smart Release 



-5- 






47 



E 45 



44 



^ 43 

o ^"^ 



*1 

41 L 



0.2 0.5 1.0 1.5 2.0 2.5 

redshift z ,, ^ 




0.0 



^ 



= pure starbursts 



A = starburst component 
of IRAS ULIRGs 



0.2 



0.4 



0.6 



log 



Fig. 3. — Starbursts with IRS measures of PAH luminosity, defined as z/L,y(7.7yum) for the 
peak f,^(7.7/im), in units of ergs s~^. Line is the formal fit to evolution of the most luminous 
sources over all redshifts, which has form (1+z)^-^. Using our adopted calibration, log Lboi = 
log z^L,^(7.7^m) - 32.80 for Lboi in L0 and z/Lj,(7.7/im) in ergs s~^; this yields that log(SFR) 
= logi/L^(7.7^m) - 42.57 for SFR in M© yr-b 

2. Evolution of the Most Luminous Starbursts and AGN 



The spectral features at rest frame wavelengths ~ 8 /im which we use to make uni- 
form measures of source luminosity can be traced with the IRS continuously with red- 
shift for < z < 3. In Figures 3 and 4, we show the results for 571 starbursts and 167 
AGN measured with the IRS. Results for starbursts are taken from the measurements and 
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Farrah et al.l (120071). Ilmanishi et al.l (12007). and ISargsvan et al. 



surements from ISargsvan fc Weedmad (120091) and ISargsvan et a 
published results from iDesai et al.l (l2009h . iDasyra et al.l ( 120091 ) 



I2OO8I): additional mea- 



(120101): an d addi tional 



Haung et al.l (120091 ). and 



Menendez-Delmestre et al.l ( 120091 ). Results for AGN are taken from the measurements and 
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Fig. 4. — All luminous AGN with dust continuum luminosity measured as z/Li,(7.9/im), in 
units of ergs s~^. Obscured AGN are those having deep silicate absorption such that the 
10 /im silicate feature absorbs more than 50% of the extrapolated unabsorbed continuum 
at 10 /xm; unobscured AGN have silicate emission. Line is the formal fit to evolution of 
the most luminous obscured AGN over all redshifts, which has form (l+z)^'®. Using our 
adopted calibration for obscured AGN, log L^oi = log z/Ly(7.9^m) - 32.63, for L^oi in Lq and 
i/Lj,(7.9/im) in ergs s~^; for unobscured AGN, the adopted calibration is log Lboi = log 
i/Lj,(7.9/im) - 32.84, indicating that bolometric luminosities are similar for obscured and 
unobscured AGN. 
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li 2003) 



j2005l ) .IWeedman et al.N2006a 



( 2008a ) 



Farrah et al. 



Weedman et al. 



. ISargsvan et al.l (|2008l). IShi et a 



(120091) 



Brand et al. 



(l2007h . JSaUna et all (l2007h. iMartinjz^Sansigreet al.l J2008f ) 



(2006b). PoUettaet al 



(12006 ) 



( 2008b ) 



(2008) 



Yan et al 



Schweitzer et al. 



( 20081). and iMarkwick-Kemper et al.l ( l2007l ): and additional published re 



suits in Hernan-Caballero et al. 



(l2009f l. 



Hao et all (l2005l ) 



The most important results in Figures 3 and 4 are that the form of luminosity evolution 
is very similar for both starbursts and AGN, and that no definitive indication is found for the 
epoch at which luminosity evolution reaches a maximum. These results mean that whatever 
process triggers the formation of luminous, dusty galaxies controls the evolution of both 
starbursts and AGN, but we cannot yet determine whether starbursts or AGN came first. 



Star Formation Rates and Obscuration of Starbursts 



Tracing the evolution of star formation in the universe is a fundamental objective of 
observational cosmology. The star formation rate (SFR) can be estimated with many ob- 



servational techniques (e.g. lKennicuttlll998l : ICalzettill2008l ). Starbursts have been observed 
to th e highest redshifts (z > 7) using rest frame ultraviolet luminosities (e.g. iBouwens et al. 



20091 ). but determination of intrinsic source luminosities and true SFRs requires large correc- 
tions for obscuration by dust. The availability of infrared spectra showing the PAH features, 
suffering much less extinction than optical or ultraviolet features, allo ws a new test of obscu- 
ration corrections (ISargsyan fc Weedmanll2009l : ISargsyan et al.ll2010l ). Results are shown in 
Figures 5 and 6. Figure 5 illustrates the dramatic contrast in SFRs that would be measured 
from the PAH 7.7 ^m feature, co r npare d to SFRs measured from the ultraviolet continuum 
(details in lSargsyan fc WeedmanI (120091 )). Even for starbursts selected in the ultraviolet, the 
differences are large. This implies large extinction for the observed ultraviolet, quantified in 
Figure 6. Using the smallest observed values of f,y(7.7 At m)/f,^(153 nm) t o dete rmine empir- 
ically the value of this ratio for unobscured starbursts, ISargsyan et al.l (120101 ) measure the 
amount by which the ultraviolet luminosity is attenuated by obscuring dust. Figure 6 shows 
that the extinction is much greater in many infrared selected starbursts, but extinction is 
still high even in ultraviolet selected starbursts. Required corrections for dust obscuration 
are higher by factors of 2 to 3 than generally assumed for the ultrav iolet selected starbursts 
used to track luminosity evolution for z > 2 (JReddy fc Steidelll2009l ) . It is crucial to confirm 
whether these large corrections for obscuration are necessary at high redshift. 



We thank Don Barry for technical assistance, Lusine Sargsyan for assistance in data anal- 
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Fig. 5. — Comparison of star formation rates measured from infrared PAH luminosity com- 
pared to observed ultraviolet luminosity; crosses are starbursts selected from Spitzer MIPS 
24 ^m surveys and triangles are starbursts selected in the ultraviolet, either as luminous 
GALEX sources or as Markarian galaxies. Luminosity z/L,^(7.7/im) is in units of ergs s~^ as 
in Figures 3 and 6. Using our adopted calibration, log L^oi = log z/L,^(7.7/im) - 32.80 for Lf,o; 
in Ijq and h'L^{7 .7 fim) in erg s^^, and log SFR(7.7/im ) = log z/L,^(7.7/im) - 42.57 for SFR in 
Mq yr~^. Ratio of deduced SFRs indicates large dust extinction of the intrinsic ultraviolet 
luminosity, shown in Figure 6. 

ysis, and V. Lebouteiller, Jeronimo Bernard-Salas and G. Sloan for developing the SMART 
optimal extraction software. Support for this work by the IRS GTO team at Cornell Uni- 
versity was provided by NASA through Contract 1257184 issued by JPL/Caltech. 
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Fig. 6. — Extinction factor for intrinsic ultraviolet luminosity, as determined from observed 
ratio fu{7-7 /im)/f,^(153 nm). Crosses are starbursts selected from Spitzer MIPS 24 /xm sur- 
veys and triangles are starbursts selected in the ultraviolet, either as GALEX sources or as 
Markarian galaxies. 
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